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Abstract 

Based on the scaling relation for the dynamics at the early time, 
a new method is proposed to measure both the static and dynamic 
critical exponents. The method is applied to the two dimensional 
Ising model. The results are in good agreement with the existing 
results. Since the measurement is carried out in the initial stage of 
the relaxation process starting from independent initial configurations, 
our method is efficient. 
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Numerical measurements of critical exponents are usually carried out from 
samples of configurations in the equilibrium generated with Monte Carlo 
algorithms. For the static exponents, Binder's method is one of the widely 
accepted [0, 0. A traditional way to obtain the dynamic exponent z is to 
measure the exponential decay of the time correlations for finite systems 
in the long-time regime 0, ||. At the critical point, except for some special 
algorithms as e.g. the cluster algorithm f| |(|, one suffers from critical slowing 
down. Recently, it has been suggested that z may be estimated from the 
power law decay of the magnetization in a relaxation process on a large 
enough lattice, in a sufficient large time but before the exponential decay 
starts 0, i- 

In the past few years, better understanding has been achieved on the 
critical relaxation processes even upto the early time of the evolution. A 
representative example of these processes is that the Ising model initially 
in random states with a small magnetization is suddenly quenched to the 
critical temperature and then evolves with the dynamics of model A. It was 
shown by Janssen et al. with e-expansion upto two-loop that besides 
the well-known long-time universal behaviour, there exists another universal 
stage of the relaxation at earlier time, termed " critical initial slip" , which 
sets in right after the microscopic time scale and eventually crosses over 
to the long-time regime. For the critical initial slip the characteristic time 
scale is to ~ tuq z ' x ° with x being a new independent critical exponent and 
m the initial magnetization. The scaling behaviour including the increase 
of the order has been illustrated by a number of authors by Monte Carlo 
simulations^ |L1], pi], [12|, |13[ , and also analytical calculations [|14], [15], [16], [17 



Of special interest is here the extension of the results in Ref. [H] to finite- 



size systems [|14| , |13[1 . In accordance to the renormalization group analysis for 
finite-size systems, we expect a scaling relation to hold for the k-th moment 
of the magnetization in the neibourhood of the critical point ]18], || |T9| . 



M {k) (t, t, L, m ) = b- k(3/u M {k) (b- z t, b 1/u r, b^L, b~ x ° 



m 



where t is the evolution time, r = (T — T c )/T c is the reduced temperature, L 
is the lattice size, and b is the spatial rescaling factor. It has been stressed in 

2 After the present letter was completed, we received a preprint of P. Grassberger 
(Wuppertal, FRG) where the critical initial slip related to damage spreading was investi- 
gated and the dynamic exponent z also measured with the Monte Carlo method in both 
two and three dimensional Ising model. 
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Ref. U that the initial states must have very short correlation lengths and 
the initial magnetization mo must be sharply prepared. The scaling relation 
for the 3-dimensional Ising model has been tested by Monte Carlo simulation 
[p~3|j . The numerical data fit into the scaling relation nicely, and Xq has been 
determined with satisfactory precision. It was observed that the microscopic 
time scale is ignorably small, for example, for the heat-bath Monte Carlo 
algorithm it is smaller than one MC sweep. This clean behaviour of the 
critical relaxation in its early time indicates a promising new way to measure 
both the static and dynamic critical exponents, which is similar to Binder's 
method in equilibrium. In the present letter, we will illustrate this idea for 
the two dimensional Ising model. 

To make the computation simpler and more efficient, we set mo to its 
fixed point, m = 0. Therefore the exponent Xq will not enter the calculation. 
Furthermore, now the time scale to — rriQ Z ^ x ° — > oo, and the critical initial 
slip gets most prominent in time direction even though the magnetization 
itself will only fluctuate around zero. The choice of mo = is essential 
in the calculation, which allows a more precise and stable measurement of 
the critical exponents. After generating randomly an initial configuration the 
system is released to the evolution with the heat-bath algorithm at the critical 
temperature. We repeat this process with independent initial configurations. 
The average is taken over the initial configurations with m = and zero 
correlation length. 

To determine z, we introduce a time- dependent Binder cumulant 

M< 4 > 

^'■'■■ £ ) 3l -3cSBS)i- (2) 

Following the scaling relation in Eq.(|l|), cumulants measured in two different 
lattices have the simple relation 

U(t,0,L 1 ) = U(b- z t,0,L 2 ) (3) 

with b = Li/L 2 . The exponent z can easily be obtained through searching 
for a time rescaling factor b~ z such that the two curves represented by both 
sides in Eq.(|3|) collapse. In other words, the cumulant can be described by a 
scaling function 

f(t/L z ) = U(t,0,L). (4) 
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Figure 1: The cumulants Ul for L — 8, 16 and 32 (labelled by +, o and 
*, respectively) together with the rescaled cumulants of double lattice (solid 
lines) . 



In the same way, it is easy to obtain other two scaling functions 

g {2 \t/L z ) = L 2l3 ^M^{t,0,L), (5) 

h(t/L z ) = L-^d T log (t, r, L) | T=0 . (6) 

With z in hand, (3/u and v can be determined by fitting respectively the 
scaling function g^(t/L z ) and h(t/L z ) from different lattices similar as in 
the determination of z. We adopt the derivative of M^ 2 '(t,r, L) for the 
estimation of u, rather than that of U(t, r, L) since the measurement of 
is more stable for large L. To our knowledge, so far the measurment of the 
static and dynamic exponents from the dynamic process in the early time has 
not been investigated. One can easily realize that at time t — > oo Binder's 
way to calculate (3jv and v is recovered^] [|TJ. Since the measurement can 
now be carried out already in the initial stage of the relaxation, the method 
becomes more efficient. 

In Fig. 1, we plot the cumulants of the two-dimensional Ising model versus 
time. The curves corresponding to L — 8, 16 and 32, repectively, are labelled 
by +,o and *. Each curve is compared with the time-rescaled cumulant 
of double lattice size given by solid lines. The exponents corresponding to 
the three pairs of best-fitted curves are z = 2.0969, 2.1493, and 2.1337, 
respectively. Our best value 2.1337 should be compared with the existing 
numerical results from J|, z = 2.13(8) and from ||20|| , z = 2.14(5) and also 



with that with the e-expansion from |]22| , z = 2.126. From the figure one can 
also see clearly that the scaling relation holds from the very beginning of the 
time evolution and how remarkably well the method is working. 

In Tab. 1, we have summarize the measured exponents z, 2(3 /v and 1/u 
for the 2-dimensional Ising model. They are in good agreement with the 



exact results [^3J and best values for z [§], |2^, p0[ . To obtain these results, 
we have averaged over 50 000 independent samples in each run and repeated 
8 runs to estimate the errors. 



3 With a similar procedure used by Binder, the critical point can also be located from 
the behaviour of U (t, r, L) around r = 0. 
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L\ <-> L 2 


2 


2/3/iz 


l/i/ 


8 «-> 16 
16 «-> 32 
32 «-> 64 


2.0969(20) 
2.1493(18) 
2.1337(41) 


0.2480(02) 
0.2494(06) 
0.2504(29) 


1.127(08) 
1.058(35) 
0.955(40) 



Table 1: Results for 2/3/z/ and 1/v, respectively, from the 2- dimensional 
Ising model. 



Compared with traditional measurements in equilibrium, much less ef- 
fort is needed with our dynamic Monte Carlo algorithm, since we do not 
enter the long-time regime where critical slowing down is severe. The ef- 
ficiency of our method may also be traced back to the simple power law 
increase of the moments at the early time. For example, it is well known 
that M( 2 )(t,0,L) ~ t(«*-W")/* when L is large enough |HJ [□]. Therefore 



the scaling relation implies M^(t, 0, L) ~ L d t^ d 2 PI V >I Z . Similar arguments 
exists for the cumulant U(t, r, L). 



In previous works on the critical initial slip H [L4], [13], [u]], the increase 
of the order and the role of m in the short- and long-time regimes are 
intensively discussed. We would like, however, to stress that even in the case 
of m = it is interesting and there exist fruitful applications, as reported 
in this letter. Our results also provide a further confirmation of the scaling 
relation discovered by Janssen et al. . Investigations for dynamics in other 
universality classes should be carried out. 
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